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Abstract 

The semi-inclusive properties of the system of neutral and charged particles with net charge equal to 
zero are considered in the grand canonical, canonical and micro-canonical ensembles as well as in micro- 
canonical ensemble with scaling volume fluctuations. Distributions of neutral particle multiplicity and 
charged particle momentum are calculated as a function of the number of charged particles. Different 
statistical ensembles lead to qualitatively different dependencies. They are being compared with the 
corresponding experimental data on multi-hadron production m. p + p interactions at high energies. 
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I. INTRODUCTION 



In relativistic high-energy colhsions many properties of produced particles follow simple 
rules of statistical mechanics. The single particle momentum spectrum approximately has the 
Boltzmann form, dN/d?p ~ exp[— (p^+m^)^/^/T], in the local rest frame of produced matter l|. 
The mean particle multiplicity of heavy particles (m ^ T) is also governed by the Boltzmann 
factor, (A^) ~ exp(— m/T). Here T, p, and m are the temperature parameter, the particle 
momentum, and its mass, respectively. The temperature parameter extracted from the data on 
p + p interactions is in the range of 160-190 MeV {2]. Thus, almost all particles are produced 
at low transverse momenta, px, and with low masses [pt, m <2 GeV). 

However, the standard statistical approach fails to reproduce the KNO scaling [Si of multi- 
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plicity distributions observed in the data on p+p, p+p, and e+ +e collisions [4, l5|, l6[ . The other 
problems are a power law behavior of the single particle transverse momentum spectrum at 
large transverse momenta, dN/d?p ~ Pj'^'', and a power law dependence of a mean^ultiplicity 
of heavy particles, (A^) ~ m"^™, where = Kp — 3 [7]. In our previous paper [8] an exten- 
sion of the standard statistical approach was suggested to a region of large transverse momenta 
and/or large particle masses [px, m > 3 GeV) by taking into account volume fluctuations. The 
proposed model, the statistical ensemble with scaling volume fluctuations (MCE/sVF), allows 
to solve the above mentioned problems of the statistical approach. 

Hadron production in high-energy collisions is characterized by two types of quantities: inclu- 
sive and semi-inclusive ones. Statistical models are usually used to describe inclusive quantities 
like mean multiplicity or mean transverse momentum. They are calculated by the summation 
over all microstates of the system with corresponding statistical weights. In the present pa- 
per we study selected properties of semi-inclusive quantities within statistical models. In this 
case the statistical summation is restricted by the additional condition, e.g. by a requirement 
that charged hadron multiplicity is fixed. The grand canonical, canonical and micro-canonical 
ensembles as well as MCE/sVF will be used. In particular, the mean multiplicity of neutral par- 
ticles and average transverse momentum of charged particles are considered at a fixed charged 
particle multiplicity. The obtained model predictions are compared to the trends observed in 
the experimental data. 

For the sake of simplicity, the system of non-interacting massless Boltzmann particles - 
neutral, positively charged and negatively charged - with the total net charge equal to zero 
Q = — N_ = is considered. The degeneracy factors are assumed to he go = g+ = g_ = 1 
and the temperature parameter is set to T = 160 MeV for quantitative calculations. 
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The paper is organized as follows. The joint multiplicity distributions of neutral and nega- 
tively charged particles, correlations, semi-inclusive averages, and the effects of quantum statis- 
tic are calculated in Section II. Semi-inclusive momentum spectra are obtained and discussed 
in Section III. A comparison with available data is presented in Section IV. Summary presented 
in Section V closes the paper. 



II. MULTIPLICITY DISTRIBUTIONS 
A. Grand Canonical Ensemble 



The grand canonical ensemble (GCE) is defined by the system volume V, temperature T, 
and charge chemical potential fig. The chemical potential fig regulates an average value of the 
conserved charge Q. For the system with zero net charge considered here fig is equal to zero. 
The mean particle multiplicities and the average energy in the GCE are: 



{No)gce = {N+)gce = {N.)g,e = N = V / 7^^ 

{E)g,, = 3T{No)gce + ST{N.),,e + 3T{N+) 



E = 9TN . 



(1) 
(2) 



In the GCE the neutral and charged multiplicities Nq and A^+, are uncorrelated and obey 
the Poisson distribution. Thus, the joint distribution of neutral A^'o and negatively charged 
particles A^_ is given by the product of two Poisson distributions which can be approximated 
by the product of two Gauss distributions at A^ ^ 1: 



PgUNo,N_) = exp{-N) X 

iVo! 



NJ. 



2N 



exp (-A^) 

X (27riV)-i/2 



(3) 
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B. Canonical Ensemble 



The canonical ensemble (CE) is described by the variables V, T, Q. The GCE expressions 
([l]l2]) for average quantities remain valid in the CE at A^ ^ 1. From the assumption Q = 
follows that A^_|_ = A^_. Consequently, the distribution of A^+ and A^_ in the CE is narrower than 
in the GCE jol. The CE distribution of neutral particles remains the same as in the GCE as it 
is not constrained by charge conservation law. The joint distribution of neutral and negatively 
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charged particles is given by 91]: 



Pce{No,N_) 



N 



exp (-A^) 



N 



■2N- 



(27riV)-i/2 



X 



exp 



/o(2iV) {NJ.y 
[No-N] 



(4) 



2N 



X (7rAr)-i/2 



exp 



(iV_ - N) 
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where Jq is the modified Bessel function. 



C. Micro- Canonical Ensemble 

The micro-canonical ensemble (MCE) is described by the variables V, E, Q. The MCE 
partition function for Nq neutral and = positively and negatively charged massless 
particles reads |lO| : 

1 1 / Y\No+2N- j^3No+6N^-l 



n 



No,N. 



(5) 



No\ (iV_)!2 Vvr2 J T{3No + 6N_) ' 

where F is the Euler gamma function. The joint probability distribution of No and in the 
MCE is, 

Qn,,n_{E,V) 



n{E,v) 



where n{E, V) = J2 



No,N- "A^o.Af- 



XE,V). 
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D. Average Multiplicities, Fluctuations and Correlations 



The mean quantities in different statistical ensembles can be expressed as 



(X) = XiNo,N.) PiNo,N.) . 

No,N- 



For the MCE distribution one obtains: 

{No)mce = (iV-)„ 



1 



{VE') 



3\l/4 



(7) 
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If the MCE energy equals the average energy of the GCE and CE, E = E, then the average MCE 
multiplicities ([H]) become equal to those in the GCE and CE ([T]). This reflects the equivalence 
of the GCE, CE, and MCE in the thermodynamic limit. However, the multiplicity distributions 
are different in these ensembles even in the thermodynamic limit. As shown in the Appendix A 
the MCE distribution dH]) for ^ 1 can be approximated as 



ttN 



exp 



{No-Nf 2{No - N){N. - N) 3{N. - Nf 



N 



N 



N 



■ (9) 
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The distributions P{No, N_) in the GCE (E]), CE (il), and MCE (ED can be written in a general 
form of the bivariate normal distribution, 

P{No,N_) ^ 



27cN -u-il- p'^) 

1 

exp 



1 f{No-Ny ^ {No-N){N.-N) , {N_ - 

2p , + 



where ui^ and u}~ are the scaled variances defined as: 



(10) 



and p is the correlation coefficient: 

0- ^ (iVoiV-) - (iVo)(iV-) ^ (iVoiV-) - (iVo)(iV_) 

V[(iV^) - m^-im - (iV_)2] . (iVo) • u;- ■ (iV_) ■ 

The averaging in Eqs. flll|12p is expressed according to Eq. ([7]). 

The scaled variances in the GCE correspond to the uncorrelated Poisson distributions ([H]): 

^Ice = ^gce = = 1 . (13) 

and the correlation coefficient fll2p is obviously equal to zero. One can similarly introduce the 
coefficients p^^ and p^~ . They are also equal to zero in the GCE. 
In the CE, 

^ce = 1 , ^ce = ^te = \ , (14) 

for the distribution (jlj). The strong correlation, = A^_, in each microscopic state of the CE 
leads to the largest possible value of the correlation coefficient: 

However, similar to the GCE, there are no correlations between neutral and charged particles, 
P'ct = 0. 

The scaled variances in the MCE are: 

Q 3 ^1 , > 

They reflect the suppression of fluctuations of neutral particles in the MCE comparing to the 
GCE and CE, and stronger suppression of fluctuations of charged particles comparing to the 
CE. The correlation coefficient, p+^g = 1 , is the same as in the CE (1151) . The exact energy 



conservation in the MCE leads to a rather strong anti-correlation between neutral and charged 
particles: 



P 



race rmce 



- ^ ^ - 0.577 



(17) 



In the large volume limit the multiplicity distribution of negatively charged particles in the 



GCE, CE, and MCE can be approximated by normal distribution 



ll|: 



P(iV_) = ^P{N_,No) = (27rcj"iV)"^/2 exp 

No 



2uj-N 



with Ljg^^ = 1, = 1/2, and cj'^g = 1/4 in the GCE, CE and MCE ensembles, respectively. 
The iV_ distribution in these ensembles is presented in Fig. [1], left. As the considered system 
has zero charge, the distributions P{N^) are equal to P{N_) ones in all statistical ensembles. 
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FIG. 1: (Color online) Examples of the multiplicity distributions (left) and the inclusive momentum 
spectra (right) of negatively charged particles obtained within the GCE, CE, MCE and MCE/sVF. 
The distributions are calculated assuming = 10 and T = 160 MeV (see text for details). 



The neutral particle multiplicity distribution reads: 



P{No) = ^P(A^_,A'o) = (27rcj°iV)~i/2 g^p 



N- 



{No -NY 
2ujm 



(19) 



with uj' 



gee = ^ee = ^ and U^^^ = 3/4. 
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E. MCE with scaling Volume Fluctuations 

The MCE with scahng volume fluctuations (MCE/sVF) [8| is described by the variables 
E, Q and V as well as the distribution function defining the scaling volume fluctuations^. All 
quantities calculated within the MCE/sVF will be denoted by the subscript a. For a description 
of the volume fluctuations it is convenient to introduce an auxiliary variable y as: 

y ^ {V/Vy/' , (20) 



and describe the scaling volume fluctuations by the scaling function ipaiv) (see Ref. 



a 



for 



details). Experimental data on the multiplicity distribution of charged hadrons in p + p inter 



actions suggest a simple analytical form of the ipaiv) function |8|, [ij, Il5| : 

My) = ^y'-'^M-ky) , (21) 

with k = 4 and T{k) being the Euler gamma function. 

The joint Nq and distribution in the MCE/sVF equals to: 

POO 

P„(iVo,iV-) = / dy P^ceiNo, N_) ^IjM , (22) 
Jo 

where PmceiNo, N_) is given by Eq. The analytical approximations for Pa{NQ,N_) are 
discussed in Appendix |Bl The inclusive mean multiplicities in the MCE/sVF are: 

{N^)a = Yl N_P^{No,N_) - AT, {No)a = Yl NoPa{No,N^) - N , (23) 

N- ,No N- ,No 

and, thus, they coincide with those in the GCE, CE, and MCE at A^ ^ 1. The inclusive 
multiplicity distributions in the MCE/sVF are: 

P^{N_) = J2MNo,N^) = = ^J^) , (24) 

No V / 

PaiNo) = J]P.(Aro,Ar_) - = (=) • (25) 

The Pq,(A^_) distribution is shown in Fig. [H It is significantly broader then the corresponding 
distributions for the GCE, CE, and MCE. The scaled variance for negatively charged and 
neutral particles is: 

= + <ce , ^- = 1^ + <ee • (26) 



^ Statistical ensembles with fluctuating extensive quantities are discussed in recent papers [l^ . 



Thus, in the MCE/sVF, because of the scahng volume fluctuations, the scaled variance increases 
in proportion to the mean multiplicity, while the scaled variance is approximately independent 
of mean multiplicity, uj ~ const, in the GCE, CE, and MCE. 

For illustration of the previously discussed properties the joint Nq and distributions 
calculated within the GCE ®, CE (HD, MCE ©, and MCE/sVF ^ are shown in Fig.H The 




N_ N 




N_ N_ 

FIG. 2: (Color online) Examples of the joint Nq and distributions calculated within the GCE 
[top left), CE (top right), MCE {bottom left) and MCE/sVF (bottom right). The distributions are 
calculated assuming = 10 (see text for details). 

multiplicities of neutral and negatively charged particles are uncorrelated in the GCE and CE 
(see Fig. [2] top panels). They are anti-correlated and correlated in the MCE and MCE/sVF, 
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respectively. A positive correlation between A^o and in the MCE/sVF is caused by the 
scaling volume fluctuations. Note that finite size effects are seen in Fig. O For example, the 
Poisson distribution is significantly asymmetric for large deviations from N. 
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FIG. 3: (Color online) Examples of the dependence of the neutral particle mean multiplicity on the 
multiplicity of negatively charged particles calculated within the GCE (top left), CE (top right), MCE 
(bottom left) and MCE/sVF (bottom right). The distributions are calculated assuming = 10 (see 
text for details). 

The distributions in Fig. [2] are rather different. Thus, it is obvious that the dependence on 

of the semi-inclusive mean multiplicity of neutral particles defined as: 

EnoNoP{No,N_) 



{Nor 



(27) 



is different in various ensembles. Namely, it is independent of in the GCE and CE: 

(iVo)* = {No):, = N. (2t 
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In the MCE, (Nq)^^^ monotonically decreases with increasing and equals approximately 
(see Appendix [C]) : 

(iV„);„ - /V (- - ^) . (29) 

Finally, the positive correlation between A^q and A^_ in the MCE/sVF leads to an approximately 
linear increase^ (No) a with increasing A^_: 

(ATo): = Ar_ . (30) 

The semi-inclusive mean multiplicities of neutral particles calculated within the GCE, CE ( 128|) . 
MCE and MCE/sVF §^ are shown as functions of A^_ in Fig. H 



F. Quantum Statistics 

In this Subsection we illustrate the effects of quantum statistics in the GCE, CE, and MCE 
using the micro-correlator method of Ref. [l6j. The mean multiplicities in the CE or MCE 
are approximately the same as in the GCE. They are given by the sum of mean occupation 



numbers with momentum p [171] 



(iV«)ee = {N'^Uce = {N^) gee - J^K) = E exp (e^/T) - 7 " ^^^^ 

where a is -|-, — , or and denotes positive, negative or neutral particles, ep = p is one particle 
energy for massless particles, 7 = +1 for Bose statistics, 7 = — 1 for Fermi statistics, and 7 = 
corresponds to the Boltzmann approximation used throughout the present paper. We study 
a neutral system, thus, chemical potentials are zero in Eq. fjHTj) . and the average occupation 
numbers (rip) = (up) are therefore the same for neutral and charged particles. The summation 
over discrete levels can be substituted by the integration in the thermodynamic limit: 

E ^ / P'dp ... . (32) 

p 







The fluctuations and correlations in the GCE, CE and MCE are very different, nevertheless they 
can be expressed in terms of the fluctuations of the occupation numbers of a single momentum 
level in the GCE, 

((A^p)').ce = v;' = {Up) (1 + 7(r^p)) • (33) 



^ A linear increase of (-/Vo)*„ce with is due to the assumption of massless particles. For non-zero value of 
mass m the relation ([50]) is changed at large N^. The maximum value of iV_ is 7V™°^ — E/2m. Eq. ((501) 
remains approximately valid for ^ jymax^ ^tVq)* approaches zero at ^max^ 
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This is a main advantage of the micro-correlator method. It allows to calculate the fluctuations 
and correlations using the following micro-correlators |18| : 



{AnlAni),,, = v^'d^^dab , (34) 



(An^Ant)ee = v;%^6a, - gV ^ '2 > (35) 



„,a2„,62 



,«2 

t^p t 

p,a p,a 



(36) 



where (5pk and 5afe are the Kronecker delta symbols, g", are particle charges, ±1 or 0, and 

1^1 - {j2<'l] ■ (E^'^") (37) 

\ p,a / \ p,a / 

is the correlation determinant. 

The variance and correlations in the GCE, CE and MCE are calculated as the sums (inte- 
grals) over momentum of the corresponding micro-correlators 



{{AN!)) = $^(An^A<) , {AN, AN,) = ^(An^A<) . (3^ 

p,k p,k 



One obtains for Bosons: 



772 



~ Ep(-P) " J-p^dp{e^/^-l)-' ~ 6C(3) - ^^'^ 



U;f^Bose ^ _ \^^^J = ^ = ^ ^± f ^ - 0.684, (40) 

^ Ep-^ _ (^p^p)' _ (^p-pS 
Ep(^p) EpW)Ep,^^g"^ EpK)Ep,-^^ 

- '^^^^U 0.407, (41) 



12C(3) 2 7r4 

where ({3) = 1.202 is the zeta Riemann function. In calculating ( l39ll4Ti) we use {n^) = {up) 
and = ^p^ = "^^p^ = for a neutral system. This gives: Epa'^p^^"^ ~ '^J2p'^p ^^"^ 
Ep a "^p^ ^p ~ 3 Ep "^p ^p- Similarly one can get the results for Fermions: 

,2 



+ ^^er™ _ ^ Q_g^2 , (42) 



9C(3) 

1 



2 



18C(3) 77r4 ^ ^ 
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The scaled variances for neutral particles are: 



, ,OBose , OBose , ,±Bose ~ -r qcq , ,OFermi , ,OFermi , ,±Fermi ~ n ni O /'/IKA 

'^gce - ^ce " ^^gce = ^■^'^^ ' ^gce " ^^ce " ^gce = U.yi/! , [4:b) 

" 603) ~ " ' "^'^^ " 903) " " • ^ ^ 

The correlation coefficients can be also calculated using micro-correlator method: 

There are no correlations between positively and negatively charged particles in the GCE, 
Pgce = O5 S'^d there is the absolute correlation in the CE and MCE, p+~ = p+~g = 1. These 
values are the same for any type of statistics. The correlation between charged and neutral 
particles, p^" = p^^ (I12p . is zero for the GCE and CE, but it has a negative value for the MCE. 
The correlation coefficient reads: 

2 



/p\ P/ A^p p p 

The Eq. ( H8l) gives for Bosons and Fermions: 



0- Base 
rmce 



-V2 



6 \ / 6 \ H -1/2 



- 2 



135C(3)2 J Vl35C(3)2 



^0— Fermi , /7) 

Pmce = - V2 



.6 



-1/2 



- 0.417 , (49) 



0.621 . (50) 



1215C(3)2 J Vl215C(3)2 

The scaled variances and correlation coeffients for the Boltzmann approximation can be 
obtained from Eqs. fl39tl41l|i8l) replacing 7 = in Eq. ( 13T1) : 

, +Boltz _ T , ,±Boltz _ l , +Boltz _ l /f-i\ 

^gce — ; ^ce ~ 2 ' ~ 4 ' 

, fiBoltz _ , fiBoltz _ 1 , fiBoltz _ n-Boltz _ ^ ~ n c;77 l'c;o^ 

'^gce - ^ce " ^ ^mce " ^ ' Pmce " ~ = - U. £)('(' . 

They, of course, coincide with our previous results fll3til4|[T"6H17l) for the Boltzmann statistics. 

One can conclude that Bose statistics makes the fluctuations always bigger and Fermi statis- 
tics - smaller: uj^^^^^ < cj^"'*^ < 00^°'^, in all statistical ensembles (GCE, CE, MCE) and for 
all types of particles (positive, negative, neutral). The strongest effect for the neutral system in 
equilibrium ^ is for the scaled variance of charged Bosons in the MCE: 'jJ^ice'^^ I ^mce^^'' — 1-6- 



^ The effects of quantum statistics for fluctuations can be much stronger at non-zero chemical potential. The 
scaled variance of Bosons may rise up to infinity near the point of Bose-condensation [lol . 
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The only correlation coefficient which fills an infiuence of quantum statistics is p in the 
MCE: p^rtce^^^^^ < Pitce^°^^^ < P^mce^"^'^ < 0- Howevcr, the quantum statistics does not change 
a sign of this correlation. Thus, the main features of the GCE, CE, and MCE fiuctuations 
and correlations found in the Boltzmann approximation - constant values of and uo^ in the 
thermodynamic limit, strong correlations, = p+~g = 1, due to the exact charge conserva- 
tion, and anti-correlation between neutral and charged particles, p[^g < 0, due to the exact 
energy conservation - remain the same for Bose and Fermi statistics. The quantum statistics 
can not simulate the MCE/sVF effects: an increase of the scaled variances in proportion to 
the mean multiplicities fl26p and a strong positive correlation, p*^^ = 1, between neutral and 
charged particles flHUl) . These new effects take place due to the scaling volume fluctuations in 
the MCE/sVF. 



III. SEMI-INCLUSIVE MOMENTUM SPECTRA 



In this section single particle momentum spectra of negatively charged particles are consid- 
ered. The inclusive spectra are denoted as F{p) and the semi-inclusive at flxed as F*{p). 
In both cases the spectra are normalized to unity: p'^dp F{p) = 1 and p'^dp F*{p) = 1. 



A. GCE and CE 



The inclusive and semi-inclusive momentum spectra in the GCE and CE are equal and read^ 
F,ce{p) = F^M = Fce{p) = F:M = ^ exp (- I) . (53) 

This follows from the fact that a single particle momentum spectrum and the particle multi- 
plicity are uncorrelated in these ensembles. 



^ This is true for the used here Boltzmann statistics. The form of momentum spectrum in the CE becomes 
different from that in the GCE for quantum gases in finite volumes. For the isospin conservation this was 



demonstrated in Ref. 



20| 
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B. Micro-Canonical Ensemble 



The inclusive single particle momentum spectrum of negatively charged particles in the MCE 
reads: 



1 1 >^ ^ AT, (3iVo + 6iV_-l)! p^No+eN.-, 
^^1^1^ (3iVo + 6iV_-4)! V- Pr.UNo,N_) , (54) 



E 



see also Ref. jiS]. The spectrum (15^ approximately has the Boltzmann form ( 153|) at momenta 
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FIG. 4: (Color online) Examples of the semi-inclusive momentum spectra of negatively charged par- 
ticles calculated within the GCE {top left), CE {top right), MCE {bottom left) and MCE/sVF {bottom 
right) for three values of N^. The distributions are calculated assuming = 10 and T = 160 MeV 
(see text for details). 



p significantly smaller than the total system energy E. However, large deviations from ([53 
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are observed close to the threshold, p = E, where the MCE spectrum approaches zero. The 
inclusive spectra F{p) in the GCE, CE (1531) . and MCE ( 15^ are shown in Fig. [H right. 

The semi-inclusive momentum spectrum at a fixed number of negatively charged particles 
is given by: 

where > 1 and C = [XIatq -Pmce(^O) ^-)] ^ is the normalization factor. Examples of the 
^mceip) spectrum for three values of are shown in Fig. H] bottom, left. The semi-inclusive 
spectra in the MCE (1551) have the Boltzmann form for p <^ E, 

KJP) = ;i exp f- , (56) 



nice \ race/ 

and the inverse slope parameter T^^^ depends on A^_. This dependence is presented in Fig.EJfor 
the MCE and other ensembles studied here. In the GCE and CE the T* is independent of 
and equal to the inverse slope parameter of the inclusive spectrum, T*^^ = T*^ = T = 160 MeV. 
In the MCE the inverse slope parameter decreases with increasing and it crosses the line 
T = 160 MeV at = A^. Thus, the inclusive momentum spectrum -Fmce(p) (15^ coincides 
with the semi- inclusive one, -F^ce(p) dSS]), at the crossing point. 

C. MCE with scaling Volume Fluctuations 

The inclusive single particle momentum spectrum in the MCE/sVF equals to: 

11^^ (3Nq + 6N_ - ly. ( px3JVo+6iv_-4 , , , 

^-M-^^EE (3A^„:6A^.-4)! (^-|) (57) 

The structure of Eq. (1371) is the same as the structure of the corresponding Eq. (13^ for the 
MCE. The only difference is in the form of the multiplicity distribution, namely Pa{NQ, A^_) is 
used in Eq. (1571) instead of Pmce{No, N_) in Eq. ( |5^ . The inclusive spectrum Fa{p) is shown 
in Fig. [1] right. It can be well approximated by the power law dependence: 

= ^'^1^^^^^^ T'^+i {p + Tk)-^-^ = 11.27 GeV^ (p + 4T)-« , (58) 

where = 4 is used in the last expression (see Ref. 

The semi-inclusive momentum spectrum in the MCE/sVF reads: 

C ~ ('3Arn + 6A^ -1)! / » x 3Afo+6Ar_-4 
^"C" - ^E j3iV:t67V:-4i! - I) ^"(^"•'^-) • <^«) 
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FIG. 5: (Color online) The dependence of the inverse slope parameter of the momentum spectra on the 
multiplicity of negatively charged particles calculated within the GCE {top left), CE (top right), 
MCE (bottom left) and MCE/sVF (bottom right). The distributions are calculated assuming = 10 
and T = 160 MeV (see text for details). 

where C = [XIatq Pa{No, N^)~\ ^. The spectrum F*(p) is plotted in Fig. Hlfor several values of 
A^_. Similar to the MCE spectrum ( l56i) the MCE/sVF one can be approximated as: 

F:{p) = ^ exp ^) , (60) 

with the inverse slope parameter T*. The dependence of T* on is shown in Fig. bottom 
right. The MCE/sVF temperature T* decreases with increasing A^_. For = the inverse 
slope parameter T* is the same in the MCE and MCE/sVF and equals to the parameter T in 
the GCE and CE. The analytical approximations of the dependence of T* on iV_ in the MCE 
and MCE/sVF are presented in Appendix O 
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IV. COMPARISON WITH DATA 



A quantitative comparison of the discussed statistical models with the experimental data 
requires a significant additional effort, which is far beyond the scope of this paper. In particular 
one should introduce proper degrees of freedom and all related conservation laws as well as a 
longitudinal collective motion of matter. Nevertheless, a qualitative comparison seems to be 
useful already, and consequently it is presented in this section. 

An excellent review of the experimental data on semi-inclusive properties of p+p interactions 



at high energies can be found in Ref. 



2l|. The volume scaling function as well as the value 



of the temperature parameter used in this work in quantitative calculations were selected in 
order to approximately reproduce the results on the inclusive distributions in p+p interactions. 
Consequently, a comparisons between these data and the model results is justified. Clearly, as 
the infiuence of global conservation laws is crucial for the considered statistical approaches the 
data referring to the semi-inclusive properties measured in the full phase-space are of primary 
importance. Several features of these data are well established {21 1. Two of them are relevant 
for the comparison with the models discussed here, namely: 

(A) the mean multiplicity of produced vr*^ mesons increases with increasing multiplicity of 
negatively charged particles and 

(B) the average transverse momentum of negatively charged particles decreases with increas- 
ing multiplicity of these particles. 

The property (B) needs several comments. First, it is well experimentally established 21| 
for charged hadron multiplicity and mean transverse momentum measured in full phase-space 
in p+p interactions at 6.6-400 GeV/c. Clearly, the full phase-space results are relevant when 
effects related to the global conservation laws are of interest. Second, the mean transverse 
momentum increases with increasing multiplicity when mid-rapidity values are considered 22 1. 
An interpretation of this dependence is, however, beyond the scope of this paper. 

For the purpose of the comparison between considered models and data the mean transverse 
momentum was calculated as: 



POO POO 

{PtY = 2 dpTpl dypTF*{p) 

Jo J -00 



—00 

00 



1 r ^2 r , / PTCOshy\ Svr^, ^ 



= dPTP^T J dypT exp(- ' j = —T* = 2.36 T* . (61) 

The mean multiplicity of neutral particles calculated within the models was identified with the 
mean tt" multiplicity. 
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GCE, CE 
MCE / sVF 



MCE 

• data 




FIG. 6: (Color online) The mean multiplicity of neutral particles (left) and the transverse momentum 
of negatively charged particles (right) as a function of the m ultip licity of negatively charged particles. 



23], left and 



24| . right are indicated by 



The experimental data on p + p interactions at 205 GeV/c 
closed circles. The predictions of the GCE, CE, MCE and MCE/sVF are shown by the lines. The 
calculations are performed assuming = 3 and T = 160 MeV (see text for details). 

The model predictions are summarized in Fig. [HI In the GCE and CE (Nq)* and (pr)* 
are independent of A^_. The MCE reproduces property (B), but leads to decrease of mean 
multiplicity (A^o)* with increasing A^_. Both features (A) and (B) are qualitatively reproduced 
in the MCE/sVF. 



V. SUMMARY 



Semi-inclusive distributions for the system of neutral and charged massless particles with 
net charge equal to zero are considered in the grand canonical, canonical, micro-canonical 
ensembles, and in the micro-canonical ensemble with scaling volume fluctuations (MCE/sVF). 
The MCE/ sVF has been included in the present study as it is the only statistical ensemble which 
reproduces the KNO scaling of multiplicity distributions and the power law behavior of the 
inclusive transverse momentum spectra measured in p + p interactions. The mean multiplicity 
of neutral particles and momentum spectra of charged particles are calculated at fixed charged 
particle multiplicity A^_. Different statistical ensembles lead to qualitatively different results 
for these semi-inclusive quantities even in the large volume limit. In other words, the semi- 
inclusive quantities can be different in different statistical ensembles despite of the ensemble 
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thermodynamical equivalence. 

The obtained model predictions are compared with the experimental data on p + p inelastic 
interactions at high energies. The MCE/sVF follows the trends observed in the data. This 
demonstrates the role of volume fluctuations in the system with exact energy and charge con- 
servation. However, the detailed comparison with the experimental data is far beyond the scope 
of this paper. The conclusive comparison with the experimental results would require inclusion 
in the statistical model calculations several neglected effects. In particular, the hadron masses 
and quantum numbers, isospin symmetry, quantum statistics, and resonance decays should be 
taken into account. 
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APPENDIX A 

The bivariate normal approximation (Q of PmceiNo, N_) can be derived as follows. Eq. ([5]) 
can be rewritten as: 

nr,„N.{E,V) = ^ exp[/(Ao,A_)], (Al) 

where 

/(Ao, A_) = (Ao + 2A_) \n[A] - ln[Ao!] - 2 ln[A_!] - ln[(3Ao + 6A_ - 1)!] , (A2) 

and A = VE^/n^. Using the Stirling formula, ln(A!) = (A + l/2)ln(A) - A + ln(27r)/2 at 
A ^ 1, the r.h.s. of Eq. ( ]A2I) can be expanded with respect to Aq and A_ near the maximum of 
/. Then, the mean multiplicities can be calculated from the condition df /ONq = df /dN_ = 0. 
Second derivatives of / with respect of Aq and A_ at the point of maximum are: 

d^f 2 d^f 6 d^f 2 

J r\_/ J rsj t/ rNJ 



dN^ A ' dNl A ' dNo ON. A 

and Eq. ([9]) follows. 
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(A3) 



APPENDIX B 



Using the approximation (Q the integration over y can be done analytically. Then the joint 
A^o and A^„ distribution in the MCE/sVF reads: 

j2 \ (fe-l)/2 

^ / -r iivoiv_ -r Oil 

r( 



dy Pmce{N^.N_) ij^iy) 



X exp{ANo + 8N_) K 



l-k 



2 J (^6 + (iV2 + 2iVoiV_ + 3m 



(Bl) 



where Ki^k is the Bessel function of the second kind. Eq. ( IBip can be simplified using the 
asymptotic expansion: 

4A;2-8A; + 3 1 



Consequently, 



TT 




2x 



e-^- 1 + 



0(x- 



X 



X > 1 . 



(B2) 



k'' 



1 fN^ + 2NoN^ + 3m 



2 \ (fc-i)/2 



X 



1 



T{k) V Qm + Nk 

Nn + 2N^ {No - AT, 



V^o + 2iV_ 



exp 



-k' 



3N 



3N_ 



(B3) 



APPENDIX C 



In the MCE with a fixed multiplicity the system temperature T^^^ can be found as 
follows. The mean multiplicity of neutral particles equals to {N^)jnce = ^^mce/^^' ^"^^ their 
average energy is (-Eo)mce = 3T^ce(^o)"ice- Thus, the total energy reads. 



^2 vj -1 • — t mce 



(CI) 



The first term in the l.h.s. of Eq. ( IClll corresponds to the average energy of neutral particles 
and the second term to that of charged particles. For the multiplicity of negatively charged 
particles close to the mean multiplicity one can solve approximately Eq. ( ICll) with respect to 
temperature. Denoting 5N_ = N_ — N and 6T = T — T^^g the solution reads: 

6N_ 



6T 



T 



2N + N 



(C2) 



or 



'4 A^_\ 



(C3) 
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Consequently, one gets: 

In the MCE/sVF at fixed one finds, 

3T: (iV*)„ + 6iV_T;, 

Using Eq. ( l30i) this gives, 

N 
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